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Using density-functional-theory calculations, we have identified new stable configurations for tri-, 
tetra-, and penta- vacancies in silicon. These new configurations consist of combinations of a ring- 
hexavacancy with three, two, or one interstitial atoms, respectively, such that all atoms remain 
fourfold. As a result, their formation energies are lower by 0.6, 1.0, and 0.6 eV, respectively, than 
the "part of a hexagonal ring" configurations, believed up to now to be the lowest-energy states. 
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Vacancies and their clusters are fundamental defects 
of silicon. Usually they result from the irradiation of 
silicon with electro ns [J B , 0, Q , neutrons d, IE 1 pro- 
tons 0, , or ions |lOlTll|. or from plastic deformations 
However, vacancy clusters can also be present 
in as-grown crystals 0] . The presence of defects in crys- 
talline semiconductors determines, to a large extent, their 
electrical and optical properties, making their study of 
great importance. 

Calculations performed using density- functional- 
theory (DFT) molecular dynamics HI E3, the 
Hartree-Fock method [H Eg, and the DFT tight- 
binding method among others, predict that the ring 
hexavacancy should be significantly more stable than 
other types of vacancy clusters. This can be explained 
using simple bond counting arguments: the crystal can 
reconstruct almost perfectly around a hexavacancy, 
making all atoms remain fourfold. For smaller clusters, 
the same calculations 0, [l^, 0] conclude that the most 
stable configurations occur when atoms are removed 
sequentially from the hexagonal ring. 

The ring hexavacancy is known to be a good trap for 
various impurities, such as carbon, oxygen, and copper 
atoms EU- It is reasonable to expect, therefore, that it 
may also be an efficient trap for self-interstitials. Explor- 
ing this avenue, we demonstrate in this Letter, on the 
basis of ab initio calculations, that penta-, tetra-, and 
tri-vacancies in the form of combinations of ring hex- 
avacancies with one, two, or three self-interstitials con- 
stitute very stable complexes, with formation energies 
significantly lower than "part of hexagonal ring" (PHR) 
configurations. In a sense, this family of defects is a 
generalization of the "fourfold coordinated point defect" 
described in ^3 , which is essentially a combination of a 
divacancy with two self-interstitials. 

The calculations of the energies and relaxed geometries 
of the vacancy clusters were performed using the Vienna 
Ab-initio Simulation Package (VASP), which employs 
pseudopotential DFT with the projector augmented- 
wave method (PAW) [HE]]. We used a 216-atom su- 
percell, an energy cutoff of 22 Ry, and the local-density 



approximation (LDA) for the exchange-correlation func- 
tional. Results are reported for F-point sampling only 
of the Brillouin zone, which we found is sufficient to en- 
sure convergence of the relative energies of the defects. 
One of the standard experimental tools for the study of 
defects in semiconductors is positron annihilation spec- 
troscopy [ilj ; we therefore also performed calculations of 
the positron lifetimes for various vacancy clusters. The 
positron wave-functions and annihilation rates were cal- 
culated using the potential and electron density given 
by the DFT calculations; the effect of electron-positron 
correlations was taken into account by introducing an ad- 
ditional correlation potential and annihilation enhance- 
ment factor according to the interpolation formulas by 
Borohski and Nieminen |22| with corrections for semi- 
conductors (see review by Puska and Nieminen |2.'"il| for 
details). 

Figure ^ presents the proposed fourfold configurations 
for the penta-, tetra-, and tri-vacancy in silicon. Figure 
n (A) (left) shows the unrelaxed ring hexavacancy, with 
labels 1 to 12 indicating the twelve atoms each having 
initially one dangling bond. In the process of relaxation, 
these twelve atoms form six new bonds with each other 
(right). However, if one silicon atom is added to the 
defect, four of the twelve atoms can form new bonds with 
it while the others pair in the same way as in the case 
of the simple hexavacancy [Fig.^(B)]. If two atoms are 
added, they will satisfy eight of the twelve dangling bonds 
while four atoms pair [Fig.Q](C)]; and if three atoms are 
added, all twelve dangling bonds of the hexavacancy are 
satisfied [Fig. U(D)]. Thus, these configurations for the 
penta-, the tetra-, and the tri-vacancy have no dangling 
bonds at all. As a result, they are expected to be more 
stable than the PHR configurations where two dangling 
bonds remain at the ends of the vacancy chain after the 
lattice has relaxed. Note that for the tetravacancy, two 
different configurations are possible: symmetric [Fig. ^ 
(C)], where the first atom is bonded to atoms 1-4 and 
the second to atoms 7 - 10, and non-symmetric, where 
the second atom is bonded to atoms 5-8. For penta- and 
tri-vacancies, all possible configurations are equivalent. 





FIG. 1: [Color online] Initial (left) and relaxed (right) geometries for the ring hexavacancy (A), and for the pentavacancy (B), 
the symmetric tetravacancy (C), and the trivacancy (D) in the fourfold configurations (combinations of a hexavacancy with 
one, two or three self-interstitials, respectively). The self-interstitials added to the hexavacancy are shown in dark green. The 
open circles indicate the positions of the atoms removed from the lattice to form the starting-point ring hexavacancy. 
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TABLE I: Calculated formation energies for various config- 
urations of vacancy clusters (N v = number of vacancies), in 
eV per defect. In the PHR configurations, atoms are removed 
sequentially from the hexagonal ring; the fourfold configura- 
tions are shown in Fig.Q For the latter, the symmetry groups 
are also given. 



N v 


Formation energy 


Group for 




PHR confs. 


Fourfold confs. 


fourfold confs. 


1 


3.51 






2 


5.01 






3 


6.80 


6.20 


D 3 


4 


8.26 (chain) 


7.26 (sym) 


C2h 




8.35 (PHR) 


7.35 (non-sym) 




5 


9.07 


8.42 


c 2 


6 


9.41 




D 3d 



The expectation that the fourfold configurations de- 
scribed above are more stable than the usual PHR con- 
figurations is verified by computing the formation ener- 
gies. In our calculations, each additional atom is initially 
placed at the geometrical center of its group of four fu- 
ture neighbours (see Fig.^J. The system is then relaxed 
using the conjugate-gradient algorithm. For comparison, 
we have also performed calculations for the PHR config- 
urations and for the chain tetravacancy. 

The calculated formation energies are presented in Ta- 
ble[J where we also give the symmetry groups of the four- 
fold defects. The formation energies for the tri-, tetra-, 
and pentavacancies in the fourfold configurations are in- 
deed lower than those for the PHR configurations, by 
0.6, 1.0, and 0.6 eV, res pect ively. It should be men- 
tioned that, in contrast to |l5j| and 0, we find the chain 
tetravacancy to have lower energy than the PHR configu- 
ration. However, both our calculations and 0,^3 show 
a very small energy difference between the two configu- 
rations so they can be considered equally stable. 

As mentioned earlier, these results were obtained using 
only the T-point to sample the Brillouin zone. In order 
to check for convergence, we have also performed some 
calculations using a 2 x 2 x 2 Monkhors-Pack grid. We find 
that the formation energies change by at most 0.5 eV, 
while the relative energies given above change by no more 
than a few percent and are thus converged with respect 
to k-point sampling. Full convergence of the formation 
energies is numerically intensive and would not alter our 
conclusions. 

Figure [3 shows the calculated binding energy, i.e., the 
energy necessary to remove one vacancy from a cluster, 
V n —* V n -\ + V. For the PHR configurations, our calcu- 
lations arc in good agreement with the results of Staab 
et al. Jjj : we also find the absolute value of the binding 
energy to be minimal for the trivacancy, and to increase 
with the size of the cluster for 3 < N v < 6. For the 




Number of vacancies N 



FIG. 2: Binding energy for vacancy clusters as a function of 
size, V n Vn-i + V . The dashed line corresponds to PHR 
configurations and the solid line to fourfold configurations (see 
Fig.0. 

fourfold configurations, now, our calculations show the 
binding energy to be approximately the same for all de- 
fect sizes. For 4 < N v < 6, this result can be explained 
by the structure of the defects: additional atoms attach 
to the hexavacancy more or less independently, and thus 
approximately the same energy is necessary to remove 
the first, second, or last atom. 

In order to identify a possible formation mechanism for 
the fourfold vacancy clusters (other than the capture of 
self-interstitials by an earlier- formed hexavacancy) , w e 
calculated, using the nudged elastic band method |24j . 
the transition barrier between the PHR and the four- 
fold pentavacancy. In this transition, an atom originally 
bonded to only two neighbours (atoms 1 and 2 in Fig. 
n (A) , for example) moves to the interstitial position to 
form two new bonds (with atoms 3 and 4 [Fig. Q] (B)]). 
Our calculations show that the barrier for this transition 
is very low, about 0.02-0.03 eV. Thus, the PHR pentava- 
cancy should quickly move to the fourfold configuration. 
In a similar way, a fourfold tetravacancy can easily be 
formed from two parallel second-nearest-neighbour diva- 
cancies by moving two twofold atoms to the interstitial 
positions. Likewise, a possible initial configuration for 
the formation of the fourfold trivacancy is three second- 
nearest-neighbour vacancies in the hexagonal ring. 

Table ILTl presents the calculated positron lifetimes for 
the different vacancy clusters. The calculations were per- 
formed for both relaxed and unrelaxed geometries. For 
the unrelaxed PHR configurations, our lifetimes are in 
perfect agreement with the results of Staab et al. [l7| . 
For the relaxed PHR states, we find close agreement for 
larger clusters, whereas for N v < 3, our calculations give 
values larger by 10-15 ps. A possible explanation for 
the (small) differences is some variations of the relaxed 
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TABLE II: Calculated positron lifetime (in ps) for the un- 
relaxed and relaxed (in brackets) geometries of the same va- 
cancy clusters as in Table [I] 



N v PHR confs. Fourfold confs. 

1 252 (226) 

2 296 (255) 

3 329 (290) 321 (258) 

4 343 (291) (chain) 342 (292) (sym) 
340 (294) (PHR) 347 (298) (non-sym) 

5 354 (301) 363 (312) 

6 376 (316) 



geometries due to the use of different methods of calcu- 
lation. 

It is generally believed that the lifetimes calculated for 
unrelaxed geometries correspond more closely to experi- 
ment because positron-induced outward relaxation com- 
pensates for the usual inward relaxation around vacancy 
clusters H3 However, in the case of the ring hex- 
avacancy, such calculations overestimate the lifetime by 
about 20 ps j2{| . Moreover, it is not clear what "unre- 
laxed geometries" means for the fourfold configurations. 
In our calculations, we simply take this as the initial con- 
figuration where the interstitial atoms are placed at the 
geometrical centers of their groups of neighbours. One 
can see from Fig. ^ that the additional atoms really 
move towards the defect center from their initial posi- 
tions in the process of relaxation, which makes our choice 
reasonable. Obviously, calculations performed for the 
geometries relaxed with respect to both electronic and 
positronic forces are necessary to get reliable values of 
the lifetimes for the fourfold configurations. Neverthe- 
less, the numbers shown in Table ITT1 provide reasonable 
estimates. 

We have performed electronic structure calculations 
and found that, like the simple ring-hexavacancy [l5lll6|. 
fourfold vacancy clusters have no energy levels in the 
band gap. As a result, they should be optically in- 
active, making their direct observation difficult. That 
being said, there is experimental evidence, mostly from 
positron annihilation spectroscopy, that the fourfold con- 
figurations are likely states of these defects. In particu- 
lar, Motoko-Kwete et al. |j| reported a positron lifetime 
value of 350 ps, consistent with both fourfold and PHR 
configurations. However, they observe the defects to be 
more stable at high temperature than the usual tetrava- 
cancies. This result has been explained by the presence 
of impurites in the material. Our calculations suggest 
that the defects actually are the fourfold configurations 
reported in this Letter. Also, the formation of fourfold 
trivacancies provides a natural explanation to the exper- 
imental results of Poirier et al. [9[; these authors have 
observed that, in the process of divacancy annealing at 



T = 250 C, the infrared absorption, which is associated 
with divacancies, decreases with time, while positron life- 
time and trapping rate remain unchanged. According 
to Table llTI the difference between positron lifetimes for 
fourfold trivacancies and divacancies is rather small, cer- 
tainly within the uncertainty arising from the computa- 
tional method (see above). Since fourfold trivacancies 
are invisible to infrared spectroscopy, the "coalescence" 
of divancies into fourfold trivacancies resolves the appar- 
ent contradiction reported in ||. 

In summary, we propose new fourfold configurations 
for tri-, tetra-, and penta- vacancies in silicon. Our DFT 
calculations show that they have formation energies lower 
by 0.6, 1.0, and 0.6 eV, respectively, than the PHR con- 
figurations, generally believed to be the stable states of 
these defects. We have identified a possible formation 
mechanism for the fourfold vacancy clusters and per- 
formed preliminary calculations for positron lifetimes as- 
sociated with them. 
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